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We calculate the single-particle Green's function of elec- 
trons that are coupled to acoustic phonons by means of higher 
dimensional bosonization. This non-perturbative method is 
not based on the assumption that the electronic system is 
a Fermi liquid. For isotropic three-dimensional phonons we 
find that the long-range part of the Coulomb interaction can- 
not destabilize the Fermi liquid state, although for strong 
electron-phonon coupling the quasi-particle residue is small. 
We also show that Luttinger liquid behavior in three dimen- 
sions can be due to quasi-one-dimensional anisotropy in the 
electronic band structure or in the phonon frequencies. 

PACS numbers: 72.10.Di, 05.30Fk, ll.10.Ef, 71.27. +a 



A famous theorem due to Migdal Q states that in cou- 
pled electron-phonon systems the electron-phonon vertex 
is for Wjj -C 1 not renormalized by phonon corrections. 
Here m is the effective mass of the electrons and M is 
the mass of the ions. An implicit assumption in the proof 
of this theorem is that the electronic system is a Fermi 
liquid |Q. However, there exists experimental evidence 
that the normal-state properties of the high tempera- 
ture superconductors do not show typical Fermi liquid 
behavior ||. Therefore it is desirable to study electron- 
phonon interactions by means of non-perturbative tech- 
niques, which do not assume a priori that the system is 
a Fermi liquid. Moreover, in heavy fermion systems the 
parameter is not necessarily small, so that Migdal's 
theorem may not be valid. In this case the self-consistent 
renormalization of the phonon energies due to the cou- 
pling to the electrons cannot be neglected j|,[3]]. 111 di- 
agrammatic approaches it is often tacitly assumed that 
the phonons remain well defined collective modes 

In this work we shall study electron-phonon interac- 
tions by means of our functional integral formulation 
of higher dimensional bosonization |7|@]j and show that 
with this non-perturbative approach one can circumvent 
the shortcomings mentioned above. Bosonization in arbi- 
trary dimensions has recently been discussed by a number 
of authors |7|-|l3|. As compared with more conventional 
operator methods Jll],[lJ] , the functional bosonization ap- 
proach j?],|]|l3) has the advantage that the retarded in- 
teraction between the electrons which is mediated via 
the phonons can be obtained trivially by integrating over 
the phonon field at the very beginning of the calcula- 
tion. Bosonization is most powerful for long-range in- 
teractions that are dominated by momentum transfers 
|q| <C kp, where kp is the Fermi momentum. In this 
case the bosonization method developed in Ref. @ leads 



to cutoff-independent results. The physically most rel- 
evant Coulomb interaction belongs to this category at 
high densities, where the Thomas-Fermi screening wave- 
vector k = (4ire 2 v) 2 is small compared with kp. Here v is 
the density of states at the Fermi energy. In this work we 
shall therefore start from a model of electrons interacting 
with Coulomb forces. For |q| <C kp the Fourier transform 
of the interaction can be approximated by its continuum 
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limit /q = -^f-- The screening problem will be solved 
explicitly by means of our bosonization approach. In this 
respect our work is complementary to studies based on 
Hubbard models JT4| , which assume that the summation 
of all diagrams describing screening processes effectively 
leads to the replacement of the long-range Coulomb in- 
teraction by a short-range Hubbard- U. 

Following the classic textbook by Fetter and Walecka 
0, we use the Debye model to describe the interac- 
tion between electrons and longitudinal acoustic (LA) 
phonons, and approximate the ionic background charge 
by a homogeneous elastic medium. Although the ions 
in real solids form a lattice, the discrete lattice struc- 
ture is unimportant for LA phonons with wave-vectors 
|q| <C kp. For a detailed description of this model and 
its physical justification see Fetter and Walecka 0). How- 
ever, some subtleties concerning screening and phonon 
energy renormalization have been ignored in Ref. H . To 
clarify these points, we first give a careful derivation of 
the effective electron-electron interaction in this model 
via functional integration. 

The dynamics of the isolated phonon system is de- 
scribed via the action 
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where b q is a complex field representing the phonons. 
Here (5 is the inverse temperature, uj m = 27tto//3, and 
q = [q, iu) m ] is a collective label. For simplicity we first 
assume LA phonons with dispersion cj q = co|q|. Here 
Co is the bare phonon velocity, which is determined by 
the short-range part of the Coulomb potential and all 
other non-universal forces between the ions. The strong 
renormalizations due to long-range Coulomb forces will 
be treated explicitly in this work. As usual, the dimen- 
sionless measure for the strength of the interaction rela- 
tive to the kinetic energy is vf^ = 



Because vf^ ^ 1 



for 



the Thomas-Fermi wave-vector k defines the 



boundary between the long- and short-range regimes. 
Note that we are implicitly assuming that the short- 
wavelength part of the problem has already been solved, 
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and that the results (for example the phonon velocity 
Co) are used as input parameters for the calculation of 
the effect of the long-wavelength modes. Such a strategy 
is also adopted in renormalization group approaches to 
condensed matter systems |l5| , and is different from the 
classic self-consistent treatment of the electron-phonon 
system by Bardeen and Pines Q. Of course, the results 
for physical quantities should be identical. We shall come 
back to this point below in our discussion of the phonon 
energy shift. 

The electronic degrees of freedom are represented by a 
Grassmann-field ip, so that the total action of the inter- 
acting electron-phonon system is 



S{ip,b} = S {Tp} + S ph {b} + S int {ip,b} . (2) 



Here 



SoW=^H^.+ek]^ (3) 



represents non-interacting spinless electrons with energy 
dispersion £k — — /i, where \x is the chemical potential. 
Cj n = 2n(n + ^)//3 is a fermionic Matsubara frequency, 
and k = [k, iCo n ]. The action Si nt {tp, b} represents the 
Coulomb energy associated with all charge fluctuations 
in the system [0 , 
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m + S^^M + SZib} , (4) 
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where V is the volume of the system, and the Fourier 
coefficients of the densities are 



Pi 



(8) 
(9) 



Here z is the valence of the ions, zN is the total number of 
conduction electrons in the system, and it is understood 
that the q = 0-term in the sums should be omitted due 
to overall charge neutrality. At this point the following 
two approximations are made in Ref. Q: (a) the bare 
Coulomb interaction in S^ ph {ip, b} is replaced by 
hand by the static screened interaction, ^§ > ^f- , and 

C[ K 

(b) the contribution Sf nt {b} is simply dropped. We shall 
see shortly that the approximation (b) amounts to ig- 
noring the self-consistent renormalization of the phonon 



frequencies |lj|5| . Although Fetter and Walecka |2| argue 
that these approximations correctly describe the physics 
of screening, it is not quite satisfactory that one has to 
rely here on words and not on calculations. Because in 
our bosonization method screening can be derived from 
first principles, we do not follow the "screening by hand" 
procedure of Ref. || , and retain at this point all terms 
in Eqs.(||)-(0) with the bare Coulomb interaction. 

We are interested in the electronic Green's function 
of the interacting many-body system. The Matsubara 
Green's function can be written as a functional integral, 



G{k) = -fl 



fV{^}V{b}e- s ^- b ^ k 4 
$V{ij)}V{b} e-BWM 
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Evidently the b- integration in Eq.(|To|) is Gaussian, so 
that it can be carried out exactly. We obtain the follow- 
ing exact expression for the interacting Green's function 



G(fc) = -p 



fV{i/)}e- s -ffW 
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The last term in Eq.(|l2|) is the effective interaction be- 
tween the electrons mediated by the phonons. Combining 
this term with S^ t {ip}, we arrive at 
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with the total effective interaction given by 
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(16) 



where the dimensionless measure for the strength of the 
electron-phonon coupling is A = ^7 2 , with 7 2 — ' 



VMu 2 ci 



ft is instructive to compare Eq.([T^) with the expres- 
sion that would result from the "screening by hand" 
procedure || described above. The approximation (a) 

amounts to setting G 2 — > in Eq. (|l2|) , while (b) 

corresponds to the replacement Wq — * 0. Then one ob- 
tains from Eq. ( |l2] ) the usual result f q — ► / q b — ■ 
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Evidently the second term can be obtained by expand- 
ing Eq. (|6|) to first order in A and replacing — > £ in 
the phonon term. By performing these replacements, one 
implicitly ignores the renormalization of the phonon en- 
ergies due to the coupling to the electrons. On the other 
hand, the effective interaction in Eq. ( |l6|) is an exact con- 
sequence of the microscopic model defined above. We 
shall see shortly that phonon energy shift and damping 
can be derived from this expression. 

Because the phonons simply modify the effective 
density-density interaction, we can obtain a non- 
perturbative expression for the Green's function by sub- 
stituting the interaction given in Eq.(|l6|) into our general 
bosonization formula for the single-particle Green's func- 
tion 0. The result for the Matsubara Green's function 
can be written as 

G(k) = V9 a (k) f dr [ dTe- i ®'- v *> r - Q " T lG a (T,T) , 
J Jo 

(17) 

where a labels "boxes" that partition the degrees of free- 
dom close Fermi surface, k a points to the center of box 
a on the Fermi surface, and the cutoff function 8 Q (k) is 
unity if k lies inside box a, and vanishes otherwise. For a 
more detailed description of this geometric construction 
see Refs. [fz|,^| Jl~0)j3~T| ] . The interacting "patch" Green's 
function G a (r~r) is of the form 



G Q (r,r) = G£(r,T)e Q ° (r ' r) 



(18) 



where Gq(t,t) is the non-interacting Green's function, 
and the Debye- Waller factor is 



Q a (r,r) 



I jRPA 1 ~ cos (q • r ~ ^t) ^ 
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Here v° = Vk£k|k=k=, and the effective interaction is 
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, U ph (q)=U a (q)+ q . 



(20) 



where IIo (q) is the usual non- interacting polarization due 
to the electronic degrees of freedom. Diagrammatically 
Eqs.(o)-(^2) are the result of a controlled resummation 
of the entire perturbation series, which is possible due to 
an underlying Ward-identity Thus, Eqs.@-@ 

contain vertex corrections to all orders in perturbation 
theory, including the renormalizations of the electron- 
phonon vertex which are only negligible if the Migdal the- 
orem is valid. 

Note that II" 1 (g) can be identified with the dressed 
phonon propagator so that it is clear that Eq.fL9|) 
takes into account that the phonon dispersion changes 
because of the coupling to the electrons. The renor- 
malized phonon mode appears as a peak in the dynamic 
structure factor ||, 



S(q,uj) = -Im ■ 

7T 1 



U ph (q,u + iO+) 

-fd>U ph (d,W + M+) 



(21) 



The qualitative behavior of S"(q, ui) can be determined 
from physical considerations In the absence of 

phonons, 5(q, u>) consists of a sum of two terms. The first 
term S co i(q, uj) is due to the collective plasmon mode. In 
d = 3 this mode approaches at long wave-lengths a finite 
value, the plasma frequency ui p i = where vp is the 
Fermi velocity. Within random-phase approximation the 



Zp'6(lu-lu p1 ) 



plasmon is not damped, so that ^^((q, uj) 

with H ZP l = f^pK^) 2 - For ^ < «i?|q| the dynamic 
structure factor has another contribution 5 sp (q, u>) due 
to the decay of density fluctuations into single-pair exci- 
tations, i.e. Landau damping. In d = 3 this is a rather 
featureless function. As long as the renormalized phonon 
velocity is small compared with vp and phonon damping 
is small, we expect that phonons give rise to an additional 
narrow peak that sticks out of the smooth background 
due to 5sp(q, w). This is the dressed phonon mode. 
Eqs.(p0[) and (21) confirm the above picture. It turns 



out that phonons are only well defined for \/X <C 
In the opposite limit the phonons are strongly damped 
and mix with the plasmon mode. For y/X <C a simple 
calculation shows that the phonon contribution to the 
dynamic structure factor can be approximated by 



S P h(q,u) 
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where the renormalized phonon energy is 
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the phonon damping is 



1 + q 2 /K 2 



A 



r = 

q 4w F |q| [l + q 2 /« 2 ] 2 ' 

and the phonon residue can be written as 
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[l+q 2 /K 2 ][l + A + q 2 /K 2 ] 



(22) 



(23) 



(24) 



(25) 



Note that at length scales small compared with the 
phonon dispersion is not renormalized, while at distances 
|q| _1 ^S> n^ 1 density fluctuations are surrounded by a 
screening cloud which modifies the phonon velocity and 
leads to phonon damping. In the long-wavelength limit 
Eq. ( p3|) implies that the phonon velocity is renormalized 
according to 



Co 



VTTx 



(26) 



For a spherical three-dimensional Fermi surface this im- 
plies in the limit of large A 



Co 
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(27) 
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This result has first been derived by Bohm and Staver 
jl7],[l8| by means of a different self-consistent treatment 
of the coupled electron-phonon system. Note that the 
renormalized phonon velocity in Eq.(p7|) is independent 
of the bare velocity cq. 

In order to investigate whether the system is a Fermi 
liquid, it is sufficient to calculate the quasi-particle 
residue Z a . The evaluation of the full momentum- and 
frequency dependent Green's function from Eqs. (]T7|)-(|2C|) 
is a formidable mathematical problem that can perhaps 
only be solved numerically. The quasi-particle residue 
associated with patch a on the Fermi surface is 0] 
Z a = e R , where R a is the constant part of the Debye- 
Waller in Eq.([f9|). Expressing f RPA in terms of the dy- 
namic structure factor and performing the Matsubara 
sum in Eq.(f9), we obtain in the limit V, f3 — > oo 



R° 



dq 

W) 1 



(/q 6 ) 2 



(uj + |v Q • q|) 2 



(28) 



We would like to emphasize that Eqs.([T7|)-(f2l|) and (J2S 
remain valid even if phonons are overdamped and mix 
with the plasmon mode. It is not difficult to see that 
the integral in Eq.(p8[) exists for arbitrary A, so that the 
Fermi liquid state is stable. To make progress analyt- 
ically, we restrict ourselves in this work to the regime 
\f\ <C where the phonon damping is negligible. Then 
we obtain R a = Rgi+Rp h , where is due to the above 
mentioned electronic terms in dynamic structure factor. 
For a spherical Fermi surface we find R" t = — if (jp) 2 : 
where r% = O(l) is a numerical constant |q]. The phonon 
contribution R^ h is obtained by substituting Eq. ( p2| ) into 

Eq. (f28[) . To leading order in VA^ the Lorentzian can 
be treated like a (5-function, so that 



R 
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Substituting the above expressions for fi q and Z? h into 
Eq. (p9[) , the integration can be performed analytically, 
with the result 



(30) 



Using (^-) = we finally obtain for the quasi- 

particle residue 



Z a = 



Vl + A 



(31) 



Note that in the limit vf —* oo (corresponding to the in- 
finite density limit) the quasi-particle residue approaches 
unity. Although our bosonization approach is only con- 
trolled for k <C hp, we expect that the result remains 
qualitatively correct at realistic metallic densities, where 

2 

-S— is of the order of unity. We conclude that at long 



wavelengths three-dimensional LA phonons always lead 
to a stable Fermi liquid, although for strong electron- 
phonon coupling the quasi-particle residue can become 
small. It should be kept in mind, however, that we have 
only retained the long-range part of the Coulomb inter- 
action, so that possible instabilities due to processes with 
large momentum transfers or superconducting pairing are 
not included in our formalism. 

It is straightforward to generalize our results for 
anisotropic systems. For example, for strictly one- 
dimensional electron dispersion the polarization in 

Eq.© is given by Ilo(g) - v ^l'JLy* - In this 
case we find that Eqs.([f7|)-(po|) give rise to Luttinger 
liquid behavior even if the phonon dispersion is three- 
dimensional. In the absence of electron-phonon interac- 
tions such a model has been studied in Ref. p9| . Al- 
ternatively, we may couple one-dimensional phonons to 
three-dimensional electrons. Then we should substitute 
in Eq.(pO|) u q = Co | | , while choosing for Tlo{q) the usual 
three-dimensional polarization. It is easy to see that in 
this case the quasi-particle residue Z a vanishes at the two 
points k" = ±kpi on the Fermi surface, where x is a unit 
vector in the x-direction. At these points |v Q -q| = vp q x , 
so that the q x -integration in Eq. ([29l) is decoupled from 
the remaining phase space. As a consequence R x is log- 
arithmically divergent. However, the total Debye- Waller 
factor Q x (r x yL, t) remains finite. (We use the label a = x 
for the patch with k" = fcpx. Because Gq (r, r) contains 
a (5-function of the components of r orthogonal to v Q , 
it is sufficient to set r — r x x and calculate Q x (r x yi, r).) 
At t = we obtain for \r x \ ^> n^ 1 



c (r x x,0) ~ -~f p h\n(K\r x \) 



with 
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The logarithmic divergence implies anomalous scaling 
characteristic for Luttinger liquids, with anomalous di- 
mension j p h- It is also easy to calculate the quasi- 
particle residue in the vicinity of the Luttinger liquid 
points ±kpii. A quantitative measure for the deviation 



For 



from these points is the parameter S = I 
^ C ^ we find that the quasi-particle residue vanishes 



|y°-*| 

|v°| 



Z a cx 



vp5 



(-0 



Iph 



(34) 



In summary, we have presented a non-perturbative ap- 
proach to the coupled electron-phonon system, which is 
not based on the assumption that the electronic sys- 
tem is a Fermi liquid. We have shown that the long- 
range Coulomb forces mediated by three-dimensional 
LA phonons cannot destabilize the Fermi liquid. How- 
ever, anisotropic electron- or phonon dispersions can 
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lead to small quasi-particle residues Z a even if the elec- 
tron dispersion is three-dimensional. Although in real- 
istic materials the phonon energies cannot be exactly 
one-dimensional on general grounds, we know ||l9|] that 
the spectral function of Fermi liquids with small quasi- 
particle residue can exhibit characteristic Luttinger liq- 
uid features in experimentally accessible regimes. More 
generally, it is tempting to speculate that the cou- 
pling between electrons and any well defined quasi-one- 
dimensional collective mode can lead to Luttinger liquid 
behavior in three-dimensional Fermi systems. In view of 
the chain-like structure of some of the high temperature 
superconductors, this result might be important for the 
explanation the unusual normal-state properties of these 
materials. 
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